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THE expansion problems which form the subject of the following 
discussion are associated with ordinary linear homogeneous differen- 
tial equations having no singular points in the interval of variation of 
the independent variable. In nearly all the problems of this sort 
that have been studied at length, the boundary conditions are of the 
type which Birkhoff! designates as regular. A notable exception 
occurs in the papers in which Liouville ? shows that a considerable 
part of the classical Sturm-Liouville theory of equations of the second 
order, including the formal expansion of a given function in series of 
characteristic functions, is capable of extension to certain special 
systems of higher order. It is readily seen that * the boundary condi- 
tions there considered are irregular. Consequently the question of 
the convergence of Liouville’s special series, which Liouville himself 
does not discuss at all, is still left open by Birkhoff’s general theory. 

The present paper is devoted to a study of this question, to which 
the author’s attention was called by Professor Bécher. No attempt 
is made to attain the greatest possible generality in the coefficients 
of the differential equation; in the matter of boundary conditions, on 
the other hand, the hypotheses are sufficiently general to include all 


1 Reference will be made to three papers of Birkhoff by means of Roman 


numerals as follows: 
: a Transactions of the American Mathematical Soeiety, 9, 219-231 
(1908). 


a wa Transactions of the American Mathematical Society, 9, 373-395 
III. Rendiconti del Circolo Matematico di Palermo, 36, 115-126 (1913). 
The definition of regular boundary conditions is introduced in II, 382-383. 
In connection with these papers of Birkhoff, see also the following articles 
by Tamarkine: 
.  Rendiconti del Circolo Matematico di Palermo, 34, 345-382 (1912). 
II. Rendiconti del Circolo Matematico di Palermo, 37, 376-378 (1914). 
21. Journal de l’Ecole Royale Polytechnique, 15, cahier 25, 85-117 (1837). 
II. Journal de Mathématiques pures et appliquées, 3, 561-614 (1838). 
8 When the differential equation is of order higher than the second. 
4 The boundary conditions are given by Liouville in non-homogeneous form, 
but it is immediately apparent that they are equivalent to a set of homogene- 
ous conditions. 
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of Liouville’s irregular cases and a large class of others. The con- 
clusion reached is not at all that suggested by the analogy of the more 
familiar expansions. It is found that functions of the sort usually 
developed in uniformly convergent series of characteristic functions 
give rise here to series which are very rapidly divergent. The explana- 
tion is that while the characteristic functions are essentially of the 
nature of real trigonometric functions in the cases usually treated, the 
corresponding functions here involve real exponential factors by which 
the order of magnitude of the terms of the series is entirely changed. 
The problem will be treated, first of all, in one of its simplest par- 
ticular cases,° which shows clearly the essential difference between the 
expansions now under consideration and the familiar ones (Section I). 
Then follows a generalization (Section IT) which is more than sufficient 
to include all of Liouville’s irregular boundary conditions. A subse- 
quent further generalization (Section III) illuminates more clearly 
the relation between these irregular cases and certain regular ones, 


I. SpectaAL DIFFERENTIAL EQUATION OF THE THIRD ORDER. 


Let us consider the differential equation ® 
(1) dx + =0 
in the interval 0 S x S 7, with the boundary conditions 


(2) u (0) = 0, u' (0) = 0, u(r) = 0. 


The parameter p is permitted to take on complex as well as real values. 
If w denotes the quantity 


5 Cf. Liouville, loc. cit. I; also H. Laudien, Dissertation, Breslau, 1914, 
pp. 67-90. In the latter memoir, with which the present writer became ac- 
quainted only after reaching the main conclusions now-expressed, there is an 
extended discussion of this special case, and it is stated that the development 
of a function satisfying the usuai conditions is uniformly convergent. Laudien 
appears to have been misled by a false analogy in thinking it allowable, in the 
demonstration of an auxiliary theorem, to dispense with a careful examina- 
tion of one of two cases, both of which are essential to the proof of the conver- 
gence of the expansion of any one arbitrary function. 

6 The discussion in the present section has numerous points of contact with 
the work of Liouville, loc. cit. I, and Laudien, loc. cit. 
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the general solution of the differential equation, when p ~ 0, has 


the form ” 
Ae— -+- Be— +. Ce — 


where A, B, and C are constants. A solution which satisfies the first 
two of the given boundary conditions for an arbitrary value of p is 
obtained by setting A= 1, B = w, C = «’, and except for an arbitrary 
constant factor this is the only solution which does so; in particular, 
not more than one characteristic function can correspond to any one 
characteristic value. A necessary and sufficient condition that 
p * 0 be a characteristic value is that the solution just mentioned 
take on the value 0 when zx = 7, that is, that p satisfy the equation 


(3) + + wie Hor = 0), 


As the problem is concerned at the outset only with the cube of p, 
there is no loss of generality in supposing that p lies in the sector of the 
complex plane for which 


—-—-<s 
(4) Sarge 


For example, the equation (3) remains unchanged if p is replaced by 
wp, the left-hand member being merely multiplied by w*. It will 
accordingly be assumed henceforth that (4) is fulfilled, unless the 
contrary is expressly stated. 

If the first term of (3) were absent, the equation would reduce to 


which can be solved explicitly by taking the logarithms of both sides. 
It has the infinitely many roots 


where 7 is any integer. Now if 9 is in the sector defined by (4), and 
large numerically, the first term of (3) is very small, while at least 
one of the others is large, and it is naturally conjectured, and readily 
proved, that throughout that part of the sector which lies outside a 
circle of sufficiently large radius about the origin the roots of (3) 


7 The general solution for p = 0 is of course A + Bx + Cz*. It is obvious 
that p = 0 is not a characteristic value, since the insertion of this expression 
in the boundary conditions gives at once A = B= C =0 


Be 
° 


386 JACKSON. 


closely correspond with those of (5). Let us show that this is the 
case. 

Let e be an arbitrarily small positive quantity, and let R be the 
region obtained by removing from the whole p-plane the interiors of all 
circles of radius e with centres at the points p’n. If the left-hand mem- 
ber of (3), minus the first term, is denoted by ¢ (p), then, in the identity 


(p) = (1 + we 


the last factor on the right is a function of p having the period 2/ v3. 
It vanishes only at the points p’n, and becomes infinite uniformly or 
approaches uniformly the limit 1 as the imaginary part of p becomes 
negatively or positively infinite. In the region R, consequently, its 
absolute value has a positive minimum 6. The absolute value of the 
parenthesis in the alternative identity 


(p) = (1 + we Vers) 


has the same minimum 6, for the two parentheses take on con- 
jugate imaginary values for conjugate imaginary values of p. Let 
a circle be described with centre at the origin and radius so large that, 
in the part of the sector (4) which lies outside it, | e-°* | < 35; and 
suppose p restricted for the moment to this part of the sector. If 
p is exterior to all the little circles, g (p) + e-°” is surely not zero, 
because of the inequalities 


(6) le | le~°r |< $6, 


the first of which is a consequence of the definition of 6 and the fact 
that at least one of the quantities e-7, e-»’e7, is in absolute value 
greater than 1. On the other hand, if p describes the circumference 
of one of the little circles in the positive sense, the argument of ¢ (p) 
increases by 27 in a single circuit, since g has just one simple ® root 
inside, and it follows from the inequalities (6) that the argument of 
¢ (p) + e-»= must do likewise. Hence (3) has just one root in each 
of the little circles. That these roots are actually and not merely 
asymptotically real, appears from the fact that complex roots of (3) 


8 The discussion of this point is put in a form adapted to the demonstration 

of Theorem I; cf. Birkhoff II, pp. 384-386; Tamarkine I, pp. 353-358. It 

ields more than is needed for the proof of Theorem II, and may be replaced 
y a simpler argument, as far as that theorem is concerned; see below. 

9 It is immediately verified that ¢() and its derivative can not vanish simul- 
taneously. The same is true of ¢(p) + e-»* for p ~ 0, so that all the roots of 
(3) which give characteristic values, even those in the finite region which has 
been temporarily left out of account, are of the first order. 
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must occur in conjugate pairs.‘° Returning to the rejected part of 
the sector (4), we need not inquire minutely what roots may be present 
there; it is sufficient for our purposes to recognize that, as ¢ is analytic 
in the finite region in question, there can not be infinitely many of 
them. There will be an index g, positive, negative, or zero, such that 
if the characteristic numbers !! are arranged in order of increasing 
absolute value, and the first is denoted by p,, the second by pgx1, 
and so on, the number pz will be nearly equal to p’n, for large values of n. 
To summarize: 

The system (1), (2) has infinitely many characteristic values on, of 
which not more than a finite number can be imaginary.!? The 
characteristic numbers can be written in the form 


1 2n 
Pn ava’ va n=9,g9+1,g9+ 
where 


To this we may add at once that there is just one characteristic 
function w(x) corresponding to each characteristic number, 


Un (x) = Pr® 1. + 
v3 


= ( VB sin — cos 


the latter expression involving only real quantities if pn is real. 

Now if the usual expansion theory were capable of extension to the 
present problem, it would be possible to develop any function, satisfy- 
ing certain not very restrictive conditions, in a uniformly convergent 
series of the form 

(8) Zan Un(2). 
But we shall see that the series is very far from possessing the necessary 
flexibility. In fact, we shall prove the following theorem: 


10 Cf. Bécher, Boundary problems in one dimension, International Congress 
of Mathematicians, Cambridge 1912; see the last paragraph of §6, containing 
a reference to Birkhoff. 

11 If some of the characteristic numbers in the finite region just mentioned 
should lie on the rays that form the boundary of the sector, we should take 
into account only the values on one of the rays. 

12 Laudien (loc. cit., pp. 73-76) gives a proof, which he ascribes to Teich- 
mann, that the characteristic values are real from the very beginning. 
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Theorem I, If a serves of the form (8) converges uniformly throughout 
any interval 


(9) 2%, 
where 0 < 2 S =, its sum must have continuous derivatives of all orders 
throughout the interval (9). 


Let x; and x2 be any two numbers subject to the inequalities 0 < 2. 
<2%1< 2%. Since the series converges uniformly throughout the 
interval (9), it converges uniformly for 


(10) 2X, 
and there must be a number G such that !° 
(11) | Qn tn (x)|<G 


for all values of x in (10) and for all values of n. 

Now if 7 is sufficiently large, pn will be real, and the interval from 
4 V3pna, to +V3pnx% will be of length greater than 27, and so will 
contain at least one value of 4 V3pax whose sine is 1. That is, there 
will be one or more values of x in (10) for which the parenthesis in (7) 
takes on the value V3. Let 2’n be such a value. Then wna(z’n) is 
positive, and 


Un (2'n) > V3 > 
Since (11) must be satisfied when 2 = 2’n, 
| dn | < Gem denn, 
On the other hand, if k is any integer, 


k 


and 
dé 
Un (a) 


< 3pn* 


for any value of x in (0, 7); we are assuming still that pn is real. If 
S m, then 


< 3 G eh — 71) 


13 Of the assumed uniform convergence, only this consequence will be used; 
and the statement of the theorem might have been made correspondingly 
more general. 
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or, in consequence of the inequalities 
< pn< 2n 


which hold for all values of n from a certain point on, 


< 3-2* Gn* 


d* 
(2) 


However large k may be, and however small numerically the negative 
quantity 22, — 2, the series which has this last expression for its 
nth term is convergent, It follows that the series obtained by differ- 
entiating (8) term by term k times is uniformly convergent * for 
0 SxS %; and as 2 may be any number between 0 and 2, the 
theorem is proved.}® 

The question is now inevitable: Is it possible to develop every 
analytic function in a uniformly convergent series of the form (8), 
with the restriction, perhaps, that the function shall vanish at the ends 
of the interval? As in other problems of this sort, it is easy to deter- 
mine what the values of the coefficients must be, if such an expansion 
exists. The formal series being given, an answer to the question may 
be expressed as follows: 

Theorem II. If N is any positive integer, however large, it 1s possible 
to define a function f(x), analytic throughout the interval OS x So, 
and vanishing with its first N derivatives both for x = 0 and for x = 7, 
such that its formal expansion in a series of the form (8) diverges at every 
untervor point of the interval. 

The system adjoint to (1), (2) consists of the differential equation 


v(0) = 0, v(r) = 0, v (r) = 0. 


Its characteristic numbers are the numbers p,, and its characteristic 
solutions have the form 


(12) = Pn(z — we wpn(xr— +. we pn(z— 


14 os course the terms of the series are continuous and in fact analytic func- 
tions of z. 

15 Professor Birkhoff points out that it may be shown by a slight modifica- 
tion of this argument that the series converges uniformly throughout a region 
of the complex z-plane including the segment (0, x2) of the axis of reals in its 
interior, and so represents a function analytic for 0 <x < Zo. 
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If a function f(x) is the sum of a uniformly convergent series (8), 
the coefficients in the series must be the numbers 


(13) an = Le 
Un(x)0n(x)dx 


The number N being given, let q be the first number greater than or 
equal to N which is congruent to 2 (mod. 3). We shall examine the 
terms of the formal series for the special function 


(14) f(x) 2(g+ 1) 


This function, a polynomial in 2, is of course analytic from 0 to 7, 
and vanishes with its first g derivatives at both ends of the interval, 
while the (q¢ + 1)th derivative takes on the value 0 at the point + 
and the value 1 at the point 0. 

To determine the order of magnitude of the numerator in (13), 
when f(x) is the function just defined, let us break up the integral 
into three, corresponding to the three terms in the right-hand member 
of (12), and consider for the moment the second of these. The index 
n may be supposed so large that pn is real and positive. Integrating 
by parts g + 2 times, we find that 


0 WPn 
— 
+ ( ) f f (g+2) (x) dz, 


WPn 


and hence, by one more integration by parts and the application of 
obvious inequalities, that 


| (a) dy — (= 


WPn 


< 
c 


prtts 


é 
where c is independent of n A corresponding inequality holds for 
the third integral, with the same value of c. As for the first, it ob- 
viously remains finite 1* as n becomes infinite, and so it is certainly 
true that 


J, (x) dz 


16 As a matter of fact, of course, it approaches zero. 
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from a certain point on. When 7 is large, therefore, 


| On (x) dx — E (=) mone | 


3c 
Pn 


e 


As q = 2 (mod. 3), the expression in brackets reduces to 


(—*) (<- en) = 92 ( *) cos = 


Pn Pn 
Now 
V3 v3 v3 
pat = cos (5 = cos (5 + sen) 


It follows that as soon as n is sufficiently large, the absolute value of 
the bracket is greater than 
3 
2pn q+ 


and hence that from a certain point on 


| J, ae| > 


With regard to the magnitude of the denominator in (13), it suffices 
to observe that, when pn is real and positive, 


| wn (x) | S Seto, | on (x) | S 


pnt 


and so 


| (x) vn (x) dx | < 


This inequality, with the one obtained for the numerator, shows 
that as n becomes infinite 


|an| >— 


where ¢; is positive and does not depend on n. 
From this relation it is apparent that the general coefficient does not 
approach zero rapidly enough to enable the terms of the series Danun(zx) 
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to remain finite as n becomes infinite. Let 2 be any particular value 
in the interval 0 < 2 < 2. Modifying the form of (7), we may write 


Un (xo) = +. 2 sin — 


The sine in the right-hand member may be zero or nearly zero for 
some values of n, but there are certainly infinitely many values of n 
for which this is not the case. For 


v3 v3 v3 
Pn4+1%0— — Pato — & ) = tot (ent 1 — &n) 
which ultimately becomes practically equal to 2, so that the two 
parentheses on the left-hand side can not both be nearly equal to 
integral multiples?’ of x. It follows that for infinitely many values 
of n 


sin( “Spur | > onto, 


| wn (ao) | > 


and 
Pn? +2 


| anuin(2o) | > 


where Co, ¢3, and cs are all positive and independent of n. As the right- 
hand member of the last inequality becomes infinite with n, the proof 
of Theorem II is complete. 

It may be remarked that no use has been made in this demonstra- 
tion of the fact that all the characteristic numbers are given by the 
asymptotic formula that we have employed. The proof would be 
unimpaired if there were infinitely many others. It is sufficient to 
know that there exist an infinite number of real characteristic values 
distributed in accordance with the formula. And this latter fact is 


17 More precisely, if 6 is the smaller of the quantities Xo, * — Xo, the index 
n may be made so large that 


Lo 9 
and then the difference between the two parentheses is between 35 and x — 46, 
and one or the other of the parentheses must itself differ from the nearest 
integral multiple of z by at least 46. © 


E€n+1 — En 


— 
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readily established by writing the characteristic equation in the real 
form 


e— or + 2eter sin 7) = 0, 


and observing the changes of sign of the left-hand member. 


Il DIFFERENTIAL EQUATION OF ORDER v WITH v — 1 SPECIALIZED 
BouNDARY CONDITIONS. 


The properties of the series studied in the preceding paragraphs are 
typical of a large class of expansion problems connected with differ- 
ential equations of order higher than the second. We shall consider 
now the differential equation 18 


15) tr @uteu=o, 23, 


where the coefficients p are real or complex functions of the real vari- 
able zx, continuous with their derivatives of all orders, in an interval 
which we may still take as the intervalO0 S237. The » linearly 
independent boundary conditions associated with this equation are 
to be subjected to the essential restriction that vy — 1 of them involve 
the values of the solution and its derivatives only at the point 0, while 
the remaining condition may involve both end-points or the point + 
alone.!® These conditions can be reduced by linear combination to the 
form 


18 It is understood that any or all of the derivatives of orders from 0 to »—2 
may occur in the equation; only the (vy — 1)th is assumed to be absent. As 
is well known, a more general differential equation can be reduced to this 
form by a change of variables; cf. Birkhoff II, p. 373, footnote. The differ- 
ential equation which Liouville considers in the paper II is specialized in a 
different way. 

19 Of course it amounts to the same thing, if a number of the conditions 
involve the derivatives at x in such a way that these can be eliminated from all 
but one of the conditions by linear combination. A similar remark is applic- 
able to the cases treated later. The conditions must include at least one 
which involves the point z, otherwise they are equivalent to the set 


u(0) = v’(0) =... = =90, 
and no characteristic solutions are pce In the system treated by Liouville, 


the rth condition involves only the point z, and the coefficients are subject to 
further restrictions. 
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W,(u) =u“ O)+ agyu”(0)=0, s=1,2,...,»—1, 
(16) 
W, (u) = ull) + (wr) + (0) = 0, 
j=0 j=0 


where 
>k,-1 2 0, y—-12k,20; 


and we shall assume them in this form. It would be possible to 
simplify them still further, if there were any object in doing so. The 
coefficients a,j, a;, @;, are any real or complex constants, and in 
particular some or all of them may vanish. 

The facts which we shall need concerning the form of solutions of the 
equation (15), apart from the question of any particular boundary 
conditions, may be taken from the papers of Birkhoff already referred 
to. Following his notation, we shall use the symbol [a], where a 
is any constant, to represent a function of p = re®, or of p and a second 
variable, which approaches the limit a uniformly with respect to 6, or 


with respect to 6 and the other variable, when r becomes infinite; the | 


same symbol will be used to represent as many different functions of 
this sort as may occur. Furthermore, the roots of the equation 
w” + 1 = Oare denoted by wi, we,...,w». Then it is possible to assign 
to each value of p a fundamental system of solutions of the equation 
(15), y,(a, p), = 1, 2,..., v, in such a way that 


(17) ys(x,p) = 


and at the same time 
k 
(18) Ys(x,p) = [1], k=1,2,...,»—1. 


These solutions have continuous derivatives of all orders with regard 
to z, foreach value of p. Iflis any one of the numbers 0,1,...,2v — 1, 
and p is restricted to the sector S, of the complex plane in which 


y 


 < argp 


the choice of this fundamental system may be so made that the func- 
tions y,(2, p) and their derivatives with regard to x shall be analytic 
functions of p for each value of 2, at least throughout the part of the 
sector exterior to a circle of sufficiently large radius about the origin. 


° 
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The last statement is still true if the sector S; is replaced by another 
one, bounded by lines parallel to the sides of S), and including the 
latter sector wholly in its interior.2° Cutting off a part of the new 
sector, if necessary, by a large circle about the origin, we will let 7, 
stand for the remaining portion, throughout which each function 
y, (x, p), with its derivatives, is analytic. Of course the regions 7, 
will overlap each other to some extent, and where they do, two differ- 
ent fundamental systems of solutions will have been defined; but as 
we shall not have occasion to use the solutions in different regions 7; 
simultaneously, no confusion will result. 

As the problem involves p only in the vth power, attention may be 
restricted to one vth of the p-plane, consisting, for example, of the 
sectors Sy and Sy-1; if there should be characteristic values on the 
boundary lines for which 


us 
arg p = 


the values on one of these lines are to be disregarded. We have to deal, 
then, only with the regions 7) and T_;, and a third region U which 
is finite in extent. It follows at once from the general existence 
theorem that in this finite region a fundamental system of solutions 
of the differential equation, and so the determinant whose vanishing 
gives the characteristic values, may be taken as analytic in p. The 
region can contain only a finite number of characteristic values, and 
these need not concern us further. 
Let us consider the region 7). Let 


wW=e”, W2=e”, 


while w3,..., w», denote the remaining vth roots of —1, arranged, for 
example, in cyclic order, though the convention as to the last v — 2 
subscripts is quite inessential.24_ Geometrically, the requirement is 
that when the roots are represented in the complex plane, as the 
vertices of a regular polygon of v sides inscribed in the unit circle, with 


20 Birkhoff III, p. 117. The assertion is readily justified by replacing the 
parameter p, in the argument of his paper I, by p + I, where I is a constant. 

21 Our use of subscripts here is different from Birkhoff’s; and in its next 
stages our discussion, though following his in general method, is vitally affected 
by the special nature of the boundary conditions. 
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a vertical side at the extreme right, the lower and upper ends of this 
side are denoted by w, and w,. respectively. Suppose that @ is a 
quantity in the interval 


0568 


The representation of the numbers w,e",..., wre, is obtained by 
rotating the original polygon through an angle @. It is clear that 
as @ varies over the interval, the argument of the complex quantity 
(w,—w,)e", where s is any subscript greater than 2, has a positive mini- 
mum. The latter statement remains true if 0 is allowed to range over 
an interval 


T 


where ¢ is a sufficiently small positive quantity. Now if p = re” is 
in 7, 8 does satisfy the last inequalities, at least when r is sufficiently 
large, and this has the consequence that as r becomes infinite, the real 
part of (w; — w,)e” remains uniformly away from zero, the real part 
of pa(w; — w,) becomes uniformly infinite, if a is any positive real 
constant, and e°*‘“s-”) uniformly approaches zero. Furthermore, 
the real part of wie” itself remains uniformly away from zero, and 
e-°“ approaches zero uniformly, even when multiplied by any posi- 
tive power of p. These facts will presently find application. 

The determinant A(p), by the vanishing of which the character- 
istic numbers are recognized, has for the element in its s-th row and 
t-th column the quantity W,(y,).. By substitution of (17) and (18) 
in (16), it is found that 


W, (yt) (pw) *s [1], s= 1, 2,° I, 
W, (yt) (pw,)*” (1] + (pwr)* (1], 


where a, is the last coefficient a; that is different from zero in (16); 
if every a; is zero, a, is understood to have the value zero likewise. 
The factor p**, which is common to all the elements of the s-th row, 
s=1, 2,..., »—1, does not affect the vanishing of A, and may 
be divided out; and we will divide out also a factor p*¥e”” from the 
elements of the last row, though this factor is not explicitly apparent 
in all the terms. The problem is thereby reduced to that of the van- 
ishing of the determinant 
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| wy" [1] we" [1] 
| [1] wok veer (we—wi) [1] 
| a,p* rw, 1] + [1] | 


w, (wy — wi) | 
[1] 


the exponent «x — /y in the last row of this determinant may be posi- 
tive, negative, or zero. 

It is at this point that the import of the remarks made at the close 
of the last paragraph but one becomes apparent. As p becomes infi- 
nite in the region 7, the first element in the last row of A; approaches 
uniformly the limit w;*”, and the elements from the third to the last 
uniformly approach zero. ‘The second of the two terms which make 
up the second element of the row is also uniformly infinitesimal. 

Let us expand the determinant according to the elements of the last 
row, and examine the cofactor of the first element. If we suppose 
each bracket symbol replaced for the moment by unity, the cofactor 
is a determinant whose value is surely different from zero, not because 
its rows consist of like powers of a set of distinct quantities (it is not 
assumed that all the exponents /, are sxccessive integers), but because, 
after division of the s-th row by wos, s = 1, 2,...,» — 1, an operation 
which of course does not affect the vanishing of the determinant, each 
column may be regarded as consisting of like powers of the v — 1 
distinct quantities the exponents of these 
powers in the several columns being the integers from 0 to v —2 
inclusive, in some order; the numbers w;/ we,...w,/ Ww, are even 
powers of we. If the brackets are restored, the cofactor is a function 
which uniformly approaches a limit different from zero when p becomes 
infinite. Similar reasoning, with a like conclusion, applies to the 
cofactor of the second element. Representing the limiting values of 
the cofactors of the first and second elements by 6; and 6 respectively, 
we see that A; has the form 


(19) Ai = [6; + [62 ; 


the terms in the expansion which involve the last »v — 2 elements 
and the second part of the second element have been merged here in 
the first bracket symbol. , 

If As is the expression which results when the bracket expressions 
in (19) are replaced by their limiting values, the equation A, = 0 can 
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be solved explicitly, having infinitely many roots given by the relation 


5; 
pn (W2 — Wi) = log (- + 2 nr, 


where the first term on the right is understood to refer to a particular 
determination of the logarithm, and nis any integer. It appears from 
the determinant expressions for 6; and 6. that each of these numbers is 
either the conjugate or the negative of the conjugate of the other; 
consequently they have the same absolute value, 6w;*» and é:w2*» 
likewise, and the above logarithm is a pure imaginary. Hence, noting 
that 


oe 
We — W, = -, 


we may write 


™ (x/v)’ 


where yo is a real constant. These roots are all of the first order, for 
the derivative of A, never vanishes. The passage from the roots of 
As to those of A; is effected by reasoning of the same sort as that used 
in connection with the simple case of the third order. The conclusion 
is as follows: 

The equation A = 0 has an infinite number of roots, which can be 
written in the form 


Mo + n+ En 


sin (x/v) ’ 


(20) Pn n= 9,9+1,94+2,...; 


where g is some integer, positive, negative, or zero,** and 73 


lim 
(2 1) n= Q. 
With the exception of a finite number at most, the roots are of the 
first order. 
Of course it still remains to consider the possibility of the existence 
of roots in the region 7_;.. But an argument analogous to that just 
completed shows that, at least from a certain point on, this region can 


22 Of course it would be possible to adjust uo so as to make g equal to zero. 
23 For the sake of convenience, a definition of €, Is given here which does not 
exactly correspond with that used in the previous special discussion. 
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contain roots only in the strip along the axis of reals which Ty_, 
has in common with 7%, so that a finite number of roots at most are 
to be added to those already found. By a change in the value of g, 
if necessary, the expression given above can be made to account for 
all the characteristic values. 

It is evident a priorz that if the coefficients in the differential equa- 
tion and in the boundary conditions are real, the conjugate of any 
characteristic value will also be a characteristic value, and it follows 
that all the characteristic values with the exception of a finite number 
will be actually and not merely asymptotically real. We do not 
restrict ourselves to this case, however. 

To each characteristic value from a certain point on corresponds 
just one characteristic function. For, 2 as we have already had 
occasion to observe, the determinant A, has first minors which are 
always different from zero when p is sufficiently large, and it follows 
that the same is true of A. An expression for the characteristic func- 
tion un(a) may be obtained by replacing the ¢-th element in the last 
row of A(pn) by y,(2, pn), t= 1,2,..., v; it is seen that the result of 
subjecting the expression so obtained to any one of the operations 
W,, 8 = 1, 2..., —1, may be written as a determinant in which 
two rows are alike, while the result of the operation W, is the vanishing 
determinant A(pn) itself. Inasmuch as a factor independent of x 
is immaterial, a characteristic solution may be obtained equally well 
by inserting the functions y; (2, pn) in place of the elements of the last 
row of A;. It is the latter more convenient form of solution that we 


24 As a matter of fact, the mere observation that the characteristic value 
is a root of the first order for A suffices to show that there can not be more 
than one characteristic function; ef., e. g., Goursat, Annales de la Faculté 
des Sciences de l|’Université de Toulouse, deuxiéme série, 10, 5-98 (1908); 
pp. 37-38; Birkhoff III, p. 118. A proof is as follows: The derivative of A 
with regard to p may be written as a sum of » determinants, in each of which 
v—1 rows are identical with the corresponding rows of A, while the elements 
of the remaining row are the derivatives of the corresponding elements of A. 
The cofactors of these derivatives in the several determinants are of course first 
minors of A, except as to algebraic sign. If p, is a characteristic value to 
which two linearly independent characteristic functions correspond, the 
first minors of A are all zero for p = p», and the derivative of A consequently 
vanishes. That is, there can not be two independent characteristic functions 
for p» unless p» is at least adouble root of A. This general fact will be assumed 
in the later discussion of other sets of boundary conditions. It is clear that 
similar reasoning shows that for a root of any order of multiplicity the number 
of linearly independent characteristic functions can never be greater than the 
order indicates; ef. Goursat, loc. cit., p. 44. I do not know whether the proof 
indicated here is to be found in the literature or not; it is the one which 
Professor Bécher gives in his lectures. 
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shall use. Or, rather, we shall use this form if the numbers 6; and 6s, 
which were introduced earlier and figure now in the coefficients of y; 
and ye, are conjugate imaginaries, but shall multiply it by 7 if 6; is 
the negative of the conjugate of 62. We have then in either case 


Vv 
(22) Un (x) [d;] 


t=1 


where d; and d:, in particular, are conjugate imaginary quantities 
different from zero. It is clear that if x is positive and n is large, all 
the terms after the second will be insignificant in comparison with 
either of the first two. Let 


d, ds ehither 
where and are real: Then 
(23) dy + dy = 2 ei + egs{ fy + pax sin (1/v)}. 


This identity is of course none the less valid for the circumstance that 
the numbers pn are not necessarily real. 

Suppose that a succession of constants 7° an,n = g,g + 1,9+2,..., 
is such that aaua(x) remains uniformly finite throughout an interval 
included in (0, 7), or, in symbols, 


(24) | dnutn (x) | << G 


for x; << 2 < where 0 << < aS 7. It is to be shown that 
throughout the interval 0 S x < a» the series 


D 
(25) antin (2) 
must converge and represent a function possessing derivatives of all 


orders. 
Because of the relations (20) and (21), the difference between the 
last factor in (23) and the expression 


(26) cos {hz + (uo +n) x} 


25 This notation is adequate only if velo the 
very beginning are simple roots of the determinant equation. In order not to 
raise questions irrelevant to the main issue, it will be assumed, in the contrary 
case, that the series (25) begins with such a value of the index n that only 
pn roots are involved. In the extension of Theorem II, where the formal 

xpansion of a definite function is under discussion, it is ‘possible to assign 
alues to the first terms of the expansion immediately; ; ef. Birkhoff II, p. 380; 
but we shall not be in any way directly concerned with these terms. 
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When n is sufficiently large, ho + (uo + n)x will vary by more than 27 
as x increases from 2; to 2, and there will surely be an intermediate 
value 2 = an’ for which the expression (26) is equal to 1. If this 
value is substituted for x in (23), for the successive values of n, the 
last factor remains from a certain point on numerically greater than a 
positive constant. Independently of this particular choice of 2, 
the quantity 


COS (w/v) / evo +n) cot — cot 


approaches zero, uniformly for 0 S 2 S x, when n becomes infinite. 


remains numerically greater than a positive constant beyond a certain 
point. Consequently the whole expression (23) remains for x = an’ nu- 
merically greater than a constant positive multiple of e+” 7*Cot(#/») | 
and the same is true of un(xn’). The inequality is strengthened if 2a’ 
in the exponent is replaced by x; it follows from (24) that 


(27) | An | < cot (a/v) 


where c’ is independent of n. 

We have now to consider the order of magnitude of the derivatives of 
un(x). Suppose x restricted to the interval 0 S x S x, where x2 < 2}. 
If k is one of the numbers 0, 1,..., »v — 1, it is at once deduced 
from (17) and (18) that 

k 
Un (2) 


< (mtn) cot (x/n) 


(28) 


where c” is independent of n. This inequality may be established 
step by step for larger values of /, by combining the inequalities 
already found with the given differential equation (15) and the equa- 
tions obtained from it by successive differentiation. Different values 
of the constant c”’ will be found for different values of k, but this is 


immaterial. 
Because of (27) and (28), the series 


d* 
an Un (a) 


converges uniformly for 0 S x S 2, and consequently represents a 
continuous function, inasmuch as the individual terms are continu- 
ous. The statement of the hypothesis for any particular value of 2; 
implies its validity for any larger value of 2; which still remains less 
than 2, and hence x2 may be any positive number less than a. With 
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the reservation that only characteristic functions corresponding to 
simple roots of the determinant equation are to be taken into account, 
we are entitled to make the assertion: 

The statement of Theorem I is applicable to the series considered in the 
present section. 

There is a corresponding extension of Theorem II. Let the left- 
hand member of the differential equation (15), without the term p’u, 
be denoted by L(u). The differential = adjoint to this is 


M (v) = (— —1 (— 


“+ p, 2, 


and the equation adjoint to (15) is 
(29) M (v) + = 0 


The determination of the adjoint boundary conditions 
carried through in some detail. 
In Lagrange’s identity, 


eL(u) — uM(v) = Plu, v), 


26 must be 


the expression P(u, v) is a bilinear form in u(x), u/(x),..., uw (2), 

v(x), v'(x),..., 0% (x). If this identity is integrated with regard to x 

from 0 to z, the resulting expression on the right, which may be de- 
zt=T 

noted by P(u, v) , 1s a bilinear form in the 27 pairs of variables 


= u(0), wm=wu'(0), u=ur-) (0), 
Up =U(T), (x), 
(0), v (0), Y= (0), 
= (7), = U(r), V2, = (zr), 


Let 7» denote the matrix of its coefficients, the coefficient of 1,2, 
standing in the s-th row and the ¢-th column. It is clear from the 
way in which 7 was obtained that the last v elements in each of 
the first v rows vanish, and likewise the first v elements in each of the 
last v rows.?” Let the linear forms W,(u),..., W.(u), which occur 


26 Cf. Birkhoff II, p. 375; Bécher, Transactions of the American Mathe- 
matical Society, 14, 403-420 (1913); pp. 404-405. Also, Liouville IT, p. 604. 

27 As a matter of fact, nearly half of the remaining elements also vanish; 
but this point is not of importance for the present discussion. 
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in (16), be supplemented by »v linear forms Wy4:(u),..., Wep(u), in 
such a way that all the 2v forms are linearly independent. If 2v 
linear forms V;(v),..., V2r(v) are determined in such a way that 


2v 


the equations 78 


V, (vw) = 0, s=v+1,v+2,..., 2», 


constitute the boundary conditions adjoint to (16). We will denote 
the matrix of the forms W by a. and that of the forms V by 2. What- 
ever these matrices may be, the matrix of the bilinear form 2W,(u) V,(v) 
is a’8, where a’ is the conjugate of a, and so the determination 
of the adjoint conditions amounts to the solution of the matrix equa- 
tion a’? = n, which gives 
B = (a’)-'n. 

Of course the reciprocal of a’ is the same as the conjugate of a-. 
Let us examine the form of the matrix /. 

The choice of the v forms Wy,1,...,W2, is arbitrary, provided that 
the condition of linear independence is fulfilled. We will choose them 
in a particularly simple way. We will set Wy: equal to that one of the 
variables w,..., uy, which does not appear as leading term in any of 
the conditions W,,..., Wy1, in (16), and will let Wy..2,..., Wey, be 
equal respectively to the variables wy1,..., Wy, with the omission of 
the one which forms the leading term of W,. It is clear that if the 
variables wu are expressed in terms of the quantities W by solving the 
linear equations which give the latter in terms of the former, the first 
v of the w’s will be completely determined by the first » — 1 equations 
and the (y+ 1)th, and so will not involve W, nor Wy,4».,..., Wo. 
That is, since the matrix of the new equations is a~, the last v — 1 
columns of a~ (and the vth column) contain no elements different 
from zero in the first y rows. The last v — 1 rows of the conjugate 
of a! contain no elements different from zero in the first v columns. 
Consequently, since the first vy columns of 7 contain non-vanishing 
elements only in the first v rows, the last v — 1 rows of the product 
(a’)—y will have all their elements in the first vy columns equal to zero. 
This means that v — 1 of the adjoint boundary conditions V = 0 


28 The use of subscripts here is different from Birkhoff’s; the change is 
made for the sake of convenience in using the language of matrices. 
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involve only the point 7 and not the point 0. The latter point 
surely will be involved in the remaining condition of the adjoint set. 
For ?° the vth row of the product matrix, like the last vy — 1 rows, 
involves only zeros in the first » columns, and it follows that the 
remaining v rows, the (v + 1)th in particular, must contain elements 
different from zero in these columns, since the determinant of the 
matrix can not vanish. The adjoint conditions may be simplified 
by linear combination to a form resembling that of the original condi- 
tions in (16). No confusion will be caused if we denote the simplified 
conditions by V,(v) = 0, s = 1, 2,..., v; the conflicting notation 
hitherto used will not be referred to again. We may write 


ki’—1 
Vi (v) = 0) + aj’ (0) + o™ (x) = 0, 
(30) 
V,(v) = + By =0, =2,3,...,%, 
where 


20, >ks’>--->k,’ 2 0. 


For each value of p, a fundamental system of solutions of the equa- 
tion (29), z,(a, p), s = 1, 2,..., v, can be assigned in such a way that 


k 
(31) (a, p) = [1], k=0,1,...,¥—1. 


It is important for us now, however, to notice that the bracket symbol 
may be given a much more precise interpretation than was done before. 
Both the relations (31) and the relations (17), (18), remain true °° 
if [a], where a is any constant, is understood to mean an expression of 
the form 
x x x E(a, p 
p p p p 

where m is an integer that may be assigned arbitrarily in advance, the 


functions ¥,(x),7 = 1, 2,..., m —1, are continuous with their de- 
rivatives of all orders for 0 S x S 7, and E(a, p) is a function of x 


29 The statement may also be proved by observing that on the contrary 
assumption the adjoint conditions could be reduced to the form v(7) = v’(r) = 
..=v@")(r) = 0, and the adjoint system could have no characteristic 


values. 
30 See Birkhoff, I, IT. 
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and p = re” which remains finite uniformly in x and 6 when r becomes 
infinite. We shall use the symbol with this meaning from now on; 
the functions y and E will be different in different cases, of course, and 
are to be replaced on occasion by quantities independent of the vari- 
able x. The solutions described in (31), regarded as functions of 2, 
are continuous with their derivatives of all orders. Their analytic 
character as functions of p need not be further specified, for the deter- 
mination of the characteristic values does not have to be repeated 
for the system (29), (30). It follows from a general theorem *} 
that the characteristic values for this system are the same as for the 
system (15), (16), and that the number of linearly independent char- 
acteristic solutions corresponding to any one characteristic number pn 
is the same for both systems. From a certain point on, the system 
(29), (30), has just one characteristic solution vn(2) for each number 
pn, and the corresponding coefficient in the formal expansion of a 


function f(x) is 
f(x) tn(a)dx 
Un(x)0n(a)dx 


We shall show how to define analytic functions f(x) for which the 
general term dnun(x) does not remain finite. 

Except for an arbitrary factor independent of x, which will be chosen 
later, the characteristic solution vn(x) may be represented by a de- 
terminant made up as follows: The ¢-th element in the first row is 
2 (a, pn), which, by (31), has the form e-°"”* [1]. The ¢-th element 
in the s-th row, s = 2, 3,..., v, is V,(z), or 


[1]. 
On division of the s-th row by (—pn)**, s = 2, 3,..., v, and multi- 


plication of the ¢t-th column by e?"”",¢ = 1, 2, ..., v, this determinant 
is reduced to the form 


(32) i, = 


[1] 1] eee [1] 


Let the limiting values of the cofactors of the first two elements of 


— 


31 See Birkhoff II, pp. 375-376; Bécher, loc. cit., p. 407. 
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the first row be denoted by 5’; and 62. As in the corresponding case 
which was discussed earlier, it is seen that neither of these numbers 
is zero, and that each is the conjugate or the negative of the conjugate 
of the other. We shall denote by u(x) in the one case the value of 
the determinant as it stands, in the other, the determinant multi- 
plied by 7. Then we have 


(34) (x) = —2), 
t=1 


where d’; and d’z are conjugate imaginary quantities different from 
zero. 

Let us examine the order of magnitude of the denominator of the 
fraction in (32), using the expressions (22) and (34) for the character- 
istic functions. In consequence of (20), (21), there is a real number 7, 
independent of n, such that 3? 


(35) R (pnw:) < (uot) cot (a/v) + 


for t= 1 and t = 2, and a fortiori for t = 3,..., v; the constant 
pocot(z/v) might be merged with y. Hence there exist numbers D, 
D’, and ¢, all independent of n, such that 


| |< Den | | < Dien(x—2) e0t 
and 


Now let us define a function f(x) by the formula (14), denoting by q 
for the present an arbitrary positive integer, and seek information 
about the corresponding value of the integral in the numerator of (32), 
with the aid of the expression (34) for m(x). By integrating by parts, 
as at the corresponding point in the preceding section, with attention 
to the special values of f(x) and its derivatives at 0 and z, and supple- 
menting this process by the application of (20), (21), and (35), it is 
seen that 


| [i@) 2) dx au (— 
0 


where ¢;, like every quantity denoted by the letter c with a subscript 
throughout the remainder of this section, is independent of n. 


Ci 
pnx Cot (x/r) 
< 


82 By R(z) is meant the real part of z. 
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It is a difference between the present problem and the special 
one treated in the preceding section that the coefficient of e*™*—” 
in (34) is not a constant, in general, but an expression of the form 


Pn Pn Pn Pn 

The additional terms, however, are readily taken into account. The 

arbitrary number m has not yet been assigned; let it be set equal to 

q+ 3. Every one of the functions f(x) ¥,(x), 7 = 1, 2,..., m — 1, 

vanishes with its first g derivatives at the points 0 and z. Hence it 

may be shown by integrating by parts that 


C2 nx COt (x/v) 
< 


Since each integral of this sort is multiplied by the reciprocal of at least 
the first power of pn, and, of course, 


we see that 


+2 | 


A corresponding inequality holds for the integral involving the 
second term of (34), the constant c, being replaced, if necessary, by a 
new constant ¢;. For t= 3, the quantity e?*”%-” remains less in 
absolute value, as m becomes infinite, than e?"”* divided by any 
power of n, and we may write 


Hence 


C7 cot (/r) 


C4 
cot (x/r) 
< 


C 
Cot t= 2,3, ...,p. 


< 


(37) 


Let d’; = where h’; and are real, so that d’o = 
By using (20) and the relation w; = cos(z/v) — 7 sin(z/v), the first 


| | 
f(x) Id i] (x—2z) dq’, | —— Qua | 
<4 
= 
| | 
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term in the braces, without the factor pn-*-*, may be written in the 
form 


(38) ent cot (x/v) ‘gt + (poten) COt (x/v) +en) 


The factor e-""*, of course, is equal to (—1)". The third factor 
approaches a limit H;— Hi as n becomes infinite. The second term 
in the braces, without the power of pn, may be written similarly as 
en™ cot (x/¥) pn=t multiplied by a factor which approaches the limit 
H, + Hx. Consequently the whole expression in braces is equal to 
(—1)" pn-2-2 multiplied by a quantity which approaches 
2H;. It follows that the absolute value of the expression is from 
some point on greater than 


C8 enn cot 
nit2 


where ¢g is a positive constant, provided that H,; #0. 

The condition that H; be equal to zero is that the imaginary part 
of the limiting value approached by the exponent of the third factor 
in (38), namely, the expression 


V 


h' ot 


be an odd multiple of $77. There is no apparent reason why this 
condition should not be satisfied, in particular cases. But it can 
surely not be satisfied for two successive values of g. We may accord- 
ingly adopt the following rule, with the assurance that the indicated 
choice of g is always possible: 

Let N be any positive integer. Then q shall be the first integer 
greater than or equal to N, for which H; + 0. 

Assuming that g is so chosen, we infer from (37) that 


da 


from some point on, and from this inequality and (36), that 


C9 


; cot (x/v) 
nit 


C10 


An > 


where ¢y and cy are positive. 
The remainder of the proof that the series diverges need not be 


elaborated. The details are suggested partly by the special discus- 
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sion in the first section, and partly by the argument just completed. 
It is found that if xp is any fixed value in the interval 0 < ap < 7, 
the absolute value of un(2) is equal to */” multiplied by a 
quantity which may be small for particular values of n, but certainly 
does not approach zero as a limit when n becomes infinite. It is 
possible to assign a positive constant ¢y, such that 


| Un (2%) | cot (2/vr) 


for infinitely many values of n. It follows that anun(2) does not 
remain finite, and 

The statement of Theorem ITI is valid for the series discussed in the 
present section. 


II]. DrrreRENTIAL EQUATION OF ORDER v WITH FEWER THAN vp — 1 
SPECIALIZED BOUNDARY CONDITIONS. 


If two of the boundary conditions, instead of a single one, involve 
the point 7, the resulting series still show essentially the same char- 
acteristics, provided that the order of the differential equation is 
greater than 4, so that the specialized conditions remain in the major- 
ity. We retain the differential equation (15) of the preceding section, 
assuming now that v 2 5, and associate with it the following boundary 
conditions: 


W, (u) = u%s) (0) + (0) = 0, 12,..,9»—2, 
(39) v—l1 


W, (u) = + Byun) + ayu CO) =0, 
O 
s=v-—l,p, 
where 
y—- 12h 0. 


The differential equation being the same, the asymptotic expressions 
(17), (18), are still available,?* and the regions S; and 7; may be de- 
fined as before. In determining the distribution of the characteristic 
values, it is sufficient to consider those in the sectors Sp and Soy_1; 
all but a finite number of these will be in the regions 7) and Tx. If 


33 For the present we may interpret the bracket symbol as it was interpreted 
at the beginning of the preceding section, without the refinement which was 
introduced later. 
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characteristic values occur on the ray arg p = 1/v, which bounds So, 
the corresponding values on the ray arg p = — 7/v may be left out of 
account. It will be convenient now to have a definite understanding 
as to four of the subscripts of the roots w,; we shall write 

", w=e", Ww=e mae’, 
so that the first four subscripts are assigned to the four right-hand 
vertices of the v-sided regular polygon which represents the v roots. 
If the polygon is rotated in the positive direction through an angle 
between 0 and the vertices representing the numbers . ., 
we, will be those furthest to the right, in the order of decreasing 
abscissas. The order of the remaining subscripts is immaterial. 

The quantities W, (y,) have the asymptotic expressions: 


(yi) = (ow)* [1], 
(yt) = (pwr) [1] + as (pwe)*s[1], s=v-—l1,p, 


where x, is the order of the highest derivative whose value at the point 0 
actually enters into the s-th condition, for the last two values of s, 
and a, is the corresponding coefficient, or is zero if the condition does 
not involve the point 0 at all. Let the factor p *s be divided from the 
s-th row of the determinant of these quantities, s = 1, 2,..., v, and 
in the last two rows let y, be written for x, — k,. Then the condition 
for a characteristic value is expressed by the vanishing of a determi- 


nant A; of the following form: *4 
[1] [1] wp [I] 
[1] wo'v-2 [1] w,*v-2 


1] + ay we [1] 
wo, [1] 
+ a, [1] 
+ a, [1] + a, p” [1] 


a, 


84 This does not correspond precisely to the determinant A, of the preceding 
section, as no exponential factor has been divided from the last rows. 


. 
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Let p be restricted to the region 7). Suppose first, for simplicity, 
that the last two of the boundary conditions do not involve the 
point 0, so that ay; and w are zero. Then in each of the minors of 
the last two rows of the determinant, an exponential factor can be 
taken from each column, and the expansion according to these minors 
has the form 


where 6; and {jj are determinants of v — 2 rows and 2 rows respec- 
tively, with powers of the quantities w; as constituents. Of these deter- 
minants, di, 613, ¢12, and {13 are surely different form zero, having 
essentially the same form as 6; and 6 in the preceding section; it is im- 
portant to observe that when the roots w, are arranged in cyclic order, w; 
and w are adjacent, and likewise w; and w3. It is immediately recog- 
nized that if p is in the remote part of 7, all the terms of (40) excepting 
possibly that with the subscripts 1, 3, are insignificant in comparison 
with the term that has the subscripts 1, 2, and may be merged with 
the latter term, by virtue of the latitude that is allowable in the 
interpretation of the bracket symbols. We may write 


(41) A; = + + wa) 


This is on the assumption that ay_; and ay are zero. Among the terms 
which have to be added to those already taken into account, if this 
restriction is removed, the worst that can occur is e?”* multiplied 
by a power of p, with a further factor which remains finite. Such a 
term will be negligible in comparison with e? + ™)*, provided that a 
power of p is negligible in comparison with e®”*, which will be the case 
if the vertex we of the polygon representing the roots remains at the 
right of the axis of imaginaries after the polygon has been rotated 
through an angle of z/v, that is, if 2x/v < 427, v > 4. As the condi- 
tion last stated was imposed at the beginning of the section, the 
terms containing powers of p may be incorporated in the expression 
(Sefp]e?" + ™*, and the validity of (41) is general. 

Suppose the bracket symbols in (41) replaced by their limiting 
values, and the resulting expression set equal to zero. Since 

W2—- W3=—e” —e =e rr 

the equation so obtained, which obviously has only roots of the first 
order, is equivalent to the relation 


(42) Qomie ” sin log (- + 
61012 


: 

= 
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where the first term on the right is a particular value of the logarithm, 
and m is any integer. As in the case previously treated, the logarithm 
is a pure imaginary. Consider, for example, the relation between 
the determinants 6 and 6,3. The latter may be regarded as obtained 
from the former by replacing ws; in the first column by we. But it 
may also be regarded as obtained by replacing each of the roots from 
w3; to we by the root next to it cyclically in the negative direction 
(whereby w3 goes into one of the roots with a subscript greater than 4, 
and w, goes into we) and then interchanging columns. This process 
amounts to a multiplication of each row of 61. by a root of unity, and 
so to a multiplication of the whole determinant by a quantity of the 
same sort. Similar reasoning applies to ¢1 and ¢13, and it follows that 
the ratio 6)3¢13/dj¢12 has its absolute value equal to unity, and the real 
part of its logarithm equal to zero. The successive values of p in 
(42) are not real, because of the factor e-*"/”, but follow each other at 
equal intervals along a ray inclined to the axis of reals at an angle of 
m/v. The correspondence between these values and the roots of the 
equation A; = 0 is established as in the earlier cases. 

The roots of A; = 0 in the region 7',_; are similarly determined. 
There will be infinitely many of them, situated asymptotically along 
the ray arg p = —7/v. But since the original problem involves 
p only in the vth power, as has already been pointed out, the existence 
of roots p along this ray implies the existence of roots pe?**/” which 
will be in the region 7», at least from a certain point on, and so will 
be among those already found. Thus the consideration of the region 
T'x_, adds only a finite number to the list of independent character- 
istic values. Hence: 

The characteristic values of the system (15), (39), are simple roots 
of the determinant equation from a certain point on, and are expres- 


sible in the form °° 
sin (27/v) n=9,9+ 1,94 2, 


where wo is real and independent of n, g is a suitable integer, positive, 
negative, or zero, and 


35 It may be that there are infinitely many characteristic values just inside 
the lower border of S.»_:, and then the formula will give instead of these the 
equivalent values in S,; but the latter will serve as well and will be more con- 
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When n is large, a characteristic solution may be obtained by setting 
p = pn in the determinant A), and replacing the ¢-th element of the 
last row by y;(2, pn) = [1], = 1, 2,...,»; and this is except for 
a constant factor the only characteristic solution °° for p = pn. The 
numbers pnw, are situated along the positive axis of reals, the numbers 
pnw, and pnw; along the rays making angles of 27/v with the positive 
axis of reals, and so on. Hence, when the determinant is expanded 
according to the minors of its last two rows, two terms are predominant 
over the others when n is large and z is not equal to 0 or 7, namely 
the terms involving and The coefficients 
of these expressions are quantities which have already been seen not 
to vanish for large values of p. The limits approached by these 
coefficients are equal in absolute value, and if the whole determinant 
is divided by either square root of their product, they will be replaced 
by limits which are conjugate imaginaries. Dividing out also the 
factor e?"™*, we may write 


(43) Un (x) = [di] + [de] 


where d; and d; are conjugate imaginary quantities different from 
zero,*® and the limits indicated by the bracket symbols are approached 
uniformly in z fora’ <2! <a" 

This expression for wn(x) is sufficient to suggest the proof of the 
first theorem, that a series of the form Lanun(x) which converges uni- 
formly throughout an interval 0 S x < 29 must represent a function 
continuous there with its derivatives of all orders. In calculating the 
order of magnitude of the terms of the derived series it will be necessary 
to go behind the relation (43) to the complete determinant expression 


venient. If the problem is real at the start, all the characteristic numbers 
from a certain point on will be actually on the ray arg p = 2/v, and not merely 
yenee it asymptotically. For if p is a characteristic number, the conjugate 
of p is a characteristic number, and e?*? ” times the conjugate of p is a character- 
istic number, and unless arg p = 2/v the number last obtained is distinct from 
p, and there are two characteristic values where, if p is large, there ought to 
be only one. It may be remarked that in the plane of the variable \ = p’, 
the characteristic values in this case are real from a certain point on, and 
become negatively infinite. 

86 Cf. footnote in connection with the corresponding passage in the preceding 
section. It is readily shown that the solution indicated can not vanish identi- 


y. 
37 When x = 0 a number of terms are of the same order of magnitude, and 
when z = x the terms involving e?"("*+™*) and ¢s(™*+%™") become com- 


parable with those indicated above. 
88 They are not the same as the numbers d; and dz of the preceding section. 
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from which it was obtained; but the argument is so similar to the 
corresponding ones already given that its details may be omitted. 

For the proof that the development of a particular analytic function 
is divergent, the form of the adjoint system is to be taken into account. 
The adjoint differential equation is of course the same as before. The 
boundary conditions may be discussed by an appropriate modification 
of the method used in the case of v—1 specialized conditions. It is 
found that v—2 of the adjoint conditions contain the point 7 only, 
while the remaining two not only may but surely do involve the 
point 0, in such a way that the sets of terms relating to this point 
in the two conditions are linearly independent.°® The conditions 
may be simplified and arranged in a manner corresponding to our 
previous practice. The characteristic functions may be written in 
the form 


(44) tn (2) = [a's] + [ay] —2), 


where the numbers pn are those already determined, d’; and d’s are 
conjugate imaginary quantities different from zero, and the approach 
of the bracket functions to their limits is uniform throughout any 
closed interval interior to (0,7). This expression, however, is not 
sufficient jor the present purpose; for in determining the order of 
magnitude of the coefficients in the series it is precisely at the point 
x = 0, where the formula breaks down, that the value of %(*) is 
needed. At this point, terms involving and 
must be explicitly taken into account. 

The characteristic function (44) was obtained from the determinant 
corresponding to*® (33) by dividing out a factor e°"”*, independ- 
ent of x, and a constant factor depending on neither pnorz. By divid- 
ing out the same exponential factor, but a different constant factor, 
a characteristic function is obtained which may be represented ade- 
quately, in a sense presently to be explained, by the expression 


Un (x) bo 63 ke’ ’ 


We ke! ePn wir 1 


39 The latter part of this statement may be proved either by making use of 
the fact that the matrix corresponding to that which was called 8 must be 
non-singular, or by observing that if the statement is assumed to be untrue 
a direct determination of the characteristic values for the adjoint system leads 
to results inconsistent with those found for the given system. 

40 This determinant in the present case has elements of the form 
[1] + B’ [1] in its second row, where k’s, and x’. are 
integers, of which y’, may be negative, and 8’: is a constant; cf. the deter- 
minant from which (43) was obtained. 
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consisting of the principal terms in the expansion of the modified 
determinant according to the minors of its first two rows; here 6’. 
and 6’; are quantities having the same absolute value, different from 
zero, and the constant factor just referred to is assumed to have been 
determined so that they are conjugate imaginaries, and k’, is the order 
of the highest derivative whose value at the point 0 enters into the 
second of the adjoint boundary conditions. Let us change our earlier 
notation to the extent of representing the new characteristic function, 
instead of (44), by en(x). The claim that this function is sufficiently 
well represented by the function ?(2) is to be understood as follows: 
Let f(x) stand for the function (14), the integer q being left arbitrary 
for the present. If the bracket symbols are given the more specific 
interpretation explained in the latter part of the preceding section, 
and the various terms are integrated by parts, it is found that 


(45) 


where c is independent of n, like the quantities to be denoted by 
the same letter with a subscript in the remainder of this section. The 
value of /f(x)tn(x) dx differs by less than ?*/”) /na+3 from 


c 


< nats 


ent cot (27/r) 


5’; 


2 
pntt 


we-4-2 
k’e 


5’. 
(46) ———— wor 


kt 
pnit We 


in comparison with which expression the difference just mentioned and 
the quantity (45) are insignificant, unless the four terms obtained by 
expanding the two determinants nearly or quite destroy each other. It 
is to be shown that this mutual destruction will not take place, at 
least if g is suitably chosen. 

If the first determinant is divided by w.-?-?+*2, and the second 
by w:-?-?+*’2, the remaining factors will be equal in the two cases, 
except as to sign, since w)/w is the same as w/w), each being equal 
to The common factor e?°¢+?)*#/” — will vanish 
only if —g-2 is congruent to k’, (mod. v), that is, for only one of 
any v consecutive values of g. On the other hand, the expression 


(47) 


pntt2 

n 
which forms the remaining factor of (46), can be treated as the expres- 
sion in braces in (37) was treated; in any v consecutive values of ¢ 
there can not be more than one for which its absolute value fails to 
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exceed a quantity /na+2, As v= 5, by hypothesis, there 
will be at least three of any v successive values of g for which 


(2) (2) de (2a/v) 


It is readily seen that 
| f Un (X) Un (x) dx 


and the rest of the divergence proof follows closely the lines of those 
already indicated. 

It suffices now to give a brief description of the principal features 
of the problem which is presented if uw of the given boundary condi- 
tions, instead of two, involve the point 7, the number y being any 
number less than half as large as v. The subscripts of the roots w; are 
to be assigned in such a way that if 0 is between 0 and 7/v the real 
parts of the numbers we” succeed each other in decreasing order of 
algebraic magnitude. In the expansion of the determinant from 
which the characteristic values are found, the largest exponentials 
will be common to the principal terms, and so will not have an influence 
on the vanishing of the function; the occurrence of roots will be ren- 
dered possible by the balancing of e?”#" against e?”#»", when yw is 
odd, and by the balancing of e°”#" against e?”#” on the one hand 
and of e?”#=7 against e?%#" on the other, when yu is even. The 
characteristic values of sufficiently large index will be roots of the first 
order, and will occur at approximately equal intervals; for odd values 
of uw they will be distributed asymptotically along the axis of reals, 
and for even values of yu, along the ray arg p = 7/v. They will be of 
amplitude 0 or 2/v exactly, from a certain point on, if the problem 
is real at the start. The characteristic functions may be represented 
by the expression 


(48) Un (x) = [di] [de] 


where d; and d, are conjugate imaginary quantities different from zero, 
and the asymptotic representation holds uniformly throughout any 
closed interval interior to (0,7). It will be observed that the numbers 
pnw, and pnw, are situated asymptotically along rays making angles 
of yux/v with the positive axis of reals, whether u is even or odd. By 
using a somewhat more explicit form for the characteristic functions, 
it may be shown that if a series of constant multiples of them con- 
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verges uniformly, its sum must necessarily have continuous deriva- 
tives of all orders. 

Of the adjoint boundary conditions, » — uw depend on the point 
m only, while the remaining yu contain linearly independent combina- 
tions of the values of v and its derivatives at the point 0. The char- 
acteristic functions of the adjoint system have an expression 


Un (x) = [ + [d’ (72), 


the interpretation of which is similar to that of (48). Finally, the 
general coefficient in the formal development of the function (14) 
may be treated as in the previous cases, and it may be shown that for 
infinitely many positive integral values *! of qg the series diverges at 
every interior point of the interval (0, 7). 

It has been pointed out that the numbers pnw, and pnw,41, which 
figure in the formula (48), are situated along rays which become 
less and less inclined to the axis of pure imaginaries as pu increases. 
If an analogous formula were to hold in the case that v is even and 
uw. = $v, which has been excluded by our hypotheses, the corresponding 
numbers would be distributed along the axis of imaginaries itself, 
and the formula would resemble those which Birkhoff finds in the case 
of regular boundary conditions.** As a matter of fact, the boundary 
conditions do become regular when yp = $7, if the order of the deriva- 
tives at the point 0 involved in the last » conditions is not too high. 
This is precisely what happens in the case of Liouville’s conditions, 
with uw = 1, when ? is equal to 2. 


41 In the formula yon megyon. to (46), there will be a determinant contain- 
Ing W1,..., Wyu-1, Wy, and one containing wW1,..., Wyu-1, When is odd, a 
common factor is disclosed by dividing the first determinant by a power of wi, 
and the second by the like power of w2; when u is even, the second determinant 
is to be divided by a power of wi, and the first by this _— of w.. Among 
v successive values of g, there will be just » — 1 values for which the determi- 
nants vanish. Not more than one other value is excluded by the requirement 
as to the order of magnitude of the quantity corresponding to (47), and there 
remain at least v — uw of every v consecutive values of g, which are available 
for the purposes of the demonstration. 

42 In the paper II, p. 389. The analogy is not perfect; the determination 
of the characteristic values becomes less simple in the regular case. 
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